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24  ABSTRACT 

25  Melioidosis,  is  an  under-reported  infectious  disease,  caused  by  the  gram-negative 

26  bacterium  Burkholderia  pseudomallei.  Understanding  disease  pathogenesis  and 

27  susceptibility  is  crucial  for  developing  newer  diagnostic  and  therapeutic  strategies  for 

28  this  life  threatening  infection.  In  this  study,  we  aimed  to  analyze  the  gene  expression 

29  levels  of  important  cytokines  in  melioidosis  patients  and  establish  useful  correlates 

30  with  disease  biomarkers  compared  to  cases  of  sepsis  infection  caused  by  other 

31  pathogens  and  healthy  individuals.  A  Qiagen  common  human  cytokines  array, 

32  profiling  the  gene  expression  of  84  important  cytokines  by  real  time  quantitative 

33  polymerase  chain  reaction  (RT  qPCR)  was  used.  We  analyzed  26  melioidosis  cases,  5 

34  healthy  controls  and  10  cases  of  sepsis  infection  caused  by  other  pathogens.  Our 

35  results  showed  a  consistent  up  regulated  expression  of  interleukins;  IL4,  IL17A, 

36  IL23A,  IL24,  interferons;  IFNAl,  IFNBl,  Tumor  necrosis  factor  (TNF)  super  family; 

37  TNFSF4,  Transforming  growth  factor  (TGF)  superfamily;  bone  morphogenetic 

38  protein  3,6  (BMP3,  BMP6),  TGFBl, other  growth  factors;  macrophage  colony 

39  stimulating  factor  (M-CSF) ,  C-fos  induced  growth  factor  (FIGF)  and  platelet  derived 

40  growth  factor  alpha  polypeptide  (PDGFA)  in  melioidosis  patients  compared  to  other 

41  sepsis  cases,  irrespective  of  comorbidities,  duration  of  fever/clinical  symptoms  and 

42  antibiotic  treatment.  Our  findings  indicate  a  dominant  Th2  and  Thl7  type  cytokine 

43  responses,  suggesting  that  their  dysregulation  at  initial  stages  of  infection  may  play  an 

44  important  role  in  disease  pathogenesis. ILIA,  ILIB  and  IL8  were  significantly  down 

45  regulated  in  septicaemic  melioidosis  patients  compared  to  other  sepsis  cases.  These 

46  differentially  expressed  genes  may  serve  as  biomarkers  for  melioidosis  diagnosis,  as 

47  targets  for  therapeutic  intervention  and  help  us  understand  immune  response 

48  mechanisms. 
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49  IMPORTANCE 

50  Melioidosis  is  a  life  threatening  infectious  disease  caused  by  a  soil-associated  gram- 

51  negative  bacterium,  B.  pseudomallei.  Melioidosis  is  endemic  in  Southeast  Asia  and 

52  northern  Australia;  however,  the  global  distribution  of  B. pseudomallei  and  the  disease 

53  burden  of  melioidosisis  still  poorly  understood.  Melioidosis  is  difficult  to  treat  as 

54  B. pseudomallei  is  intrinsically  resistant  to  many  antibiotics  and  requires  a  long  course 

55  of  antibiotic  treatment.  Mortality  rate  remains  high  in  endemic  areas  with 

56  reoccurrence  being  common.  Therefore,  it  is  imperative  to  diagnose  the  disease  at  an 

57  early  stage  and  provide  vital  clinical  care  to  reduce  the  mortality  rate.  With 

58  limitations  in  treatment  and  lack  of  a  vaccine,  it  is  crucial  to  study  the  immune 

59  response  mechanisms  to  this  infection  to  get  a  better  understanding  of  disease 

60  pathogenesis  and  susceptibility.  Therefore,  this  study  aimed  to  analyze  the  gene 

61  expression  levels  of  important  cytokines  to  establish  useful  correlation  for  diagnostic 

62  and  therapeutic  purposes. 

63  INTRODUCTION 

64  Melioidosis  is  a  life  threatening  infectious  disease  and  is  endemic  in  Southeast  Asia 

65  and  northern  Australia  (1).  A  recent  report  estimates  melioidosis  disease  burden  to  be 

66  165,000  cases  per  year  (2).  Lack  of  awareness  of  melioidosis  disease  among 

67  physicians  and  lack  of  diagnostic  methods  contribute  to  underreporting  in  many 

68  endemic  countries.  Infection  is  suspected  to  be  acquired  mainly  via  skin  during 

69  exposure  to  soil  and  contaminated  water.  Nevertheless  inhalation  of  aerosolized 

70  bacteria  during  extreme  weather  events  such  as  rainfall  and  storms  has  also  been 

71  reported  (2,  3). The  disease  is  strongly  associated  with  comorbidities  such  as  diabetes 

72  mellitus,  chronic  kidney  disease,  thalassemia,  immunosuppression  and  excessive 
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73  alcohol  intake  (1,  4,  5).  A  broad  spectrum  of  clinical  presentations  ranging  from  acute 

74  fulminant  septicemia  to  chronic  localized  abscesses  are  reported  for  melioidosis  (5). 

75  Early  diagnosis  and  appropriate  antibiotic  treatment  plays  a  crucial  role  in  preventing 

76  mortality  and  recurrence.  Advancement  of  new  immunodiagnostic  methods  and 

77  therapeutic  strategies  is  important  for  disease  management  of  melioidosis,  given  the 

78  lack  of  vaccines  and  limitations  in  drug  treatment  (3). 

79  Studying  the  host  immune  responses  to  infection  is  crucial  for  understanding  disease 

80  pathogenesis,  susceptibility  and  immune  correlates  of  protection  (3).  Cytokines  are 

81  vital  immune  modulators  that  regulate  and  determine  the  nature  of  immune  responses 

82  to  an  infection  (6).  Activation  of  leukocytes  and  cytokine  networks  are  prominent 

83  features  of  inflammation  and  the  septic  response  (7).  Pro-  and  anti-inflammatory 

84  cytokines  play  a  critical  role  in  regulating  overall  immune  responses  and  in 

85  establishing  homeostasis,  and  their  dysregulation  is  instrumental  in  triggering  disease 

86  progression  and  severity  (8).  Hence  a  detailed  study  of  the  cytokine  cascade  events  at 

87  the  transcriptome  level  during  an  infection  is  useful  to  understand  disease 

88  pathogenesis  and  susceptibility.  Although  cytokine  cascade  events  following 

89  B. pseudomallei  infection  have  been  studied  in  several  animal  models  (7,  9-13),  data 

90  on  human  host  mRNA  expression  levels  of  cytokines  is  limited.  Pro-inflammatory 

91  cytokines  such  as  IL8,  IL6,  IL12,  IL18,  IL15,  IFNy,  TNFa,  ILip,  anti-inflammatory 

92  cytokines  such  as  IL4  and  several  other  chemokines  have  been  implicated  in  disease 

93  outcome  during  the  early  acute  phase  of  B. pseudomallei  infection  (7,  14-16).While 

94  individual  cytokines  have  been  investigated  in  previous  studies^the  profiling  of  entire 

95  cytokine  networks  is  necessary  to  comprehensively  understand  specific  immune 

96  response  pathways  and  thereby  the  pathophysiology  of  melioidosis.  Such  a  profile 

97  may  also  help  identify  disease  biomarkers  with  therapeutic  implications. 
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98  We  have  successfully  established  a  nation-wide  surveillance  system  in  Sri  Lanka 

99  which  has  resulted  in  finding  more  confirmed  cases  of  melioidosis  (unpublished  data). 

100  In  this  study  we  aimed  to  analyze  the  gene  expression  profiles  of  important  human 

101  cytokines  in  Sri  Lankan  melioidosis  patients  to  further  understand  the  immune 

102  response  mechanisms  during  melioidosis  and  establish  useful  correlates  with  disease 

103  biomarkers. 

104  METHODS 

105  Patient  enrollment 

106  Nationwide  active  surveillance  for  melioidosis  was  established  in  multiple  state  and 

107  private  hospitals  throughout  Sri  Lanka,  with  ethics  approval  from  the  Ethics  Review 

108  Committee,  Faculty  of  Medicine,  University  of  Colombo,  Sri  Lanka  and  the  Office  of 

109  Human  Research  Protection  (OHRP),  United  States  Army  Medical  Research  and 

110  Material  Command  (USAMRMC).  Patients  fitting  the  clinical  case  definition  of 

111  melioidosis  i.e.  febrile  illness  for  more  than  5  days,  pneumonia,  septic  arthritis,  skin 

112  lesions,  septicaemia,  lung,  soft  tissue  or  deep  abscess  were  recruited  for  initial 

113  screening  for  melioidosis.  Blood,  pus  and  other  patient  specimens  were  collected  for 

114  bacterial  cultures  and  serum  samples  were  collected  for  indirect  haemagglutination 

115  (IHA)  antibody  test.  Any  positive  bacterial  cultures  were  further  screened  and 

116  confirmed  as  B.pseudomallei  by  PCR.  All  samples  for  the  study  were  collected 

117  between  September  2014  and  April  2016. 

118  Patients  who  were  culture  positive  for  B.  pseudomallei  and  /  or  had  high  antibody 

119  titers  (>640)  by  the  IHA  test  were  recruited  for  our  study  and  classified  as  positive 

120  cases  of  melioidosis.  Culture  and  PCR  positive  samples  were  considered  as  confirmed 

121  cases  of  melioidosis.  Samples  with  an  antibody  titre  of  >640  by  IHA  testing  were 
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122  considered  as  probable  eases  of  melioidosis.  At  the  time  of  recruitment  all  melioidosis 

123  patients  were  undergoing  antibaeterial  treatment. 

124  We  also  reeruited  healthy  donors  and  patients  fitting  the  elinical  definition  of  severe 

125  sepsis/septic  shock  (as  per  the  2012  WHO  guidelines  for  sepsis  management)  who 

126  were  negative  for  B.  pseudomallei,  as  negative  controls  for  our  gene  expression 

127  profiling  study  (17). 

128  Bacterial  culture  and  identification 

129  Primary  isolation  of  B.  pseudomallei  was  done  at  the  admitting  hospital  using 

130  conventional  culture  techniques  for  blood,  sputum,  pus  and  other  specimens.  Bacterial 

131  isolates  that  were  oxidase  positive,  gentamicin-resistant  and  gram-negative  bacilli 

132  were  forwarded  to  the  reference  laboratory  in  Colombo  where  they  were  sub-eultured 

133  to  establish  pure  growth  and  maintained  at  -70°C  in  15%  brain  heart  infusion  (BHI) 

134  glyeerol  for  subsequent  definitive  tests.  Baeteria  were  resuscitated  by  subculture  onto 

135  5%  blood  agar  and  ineubated  for  24  h  at  37  °C  to  give  single  eolony  growth  for  all 

136  subsequent  tests. 

137  Real  time  PCR  assay  for  confirmation  of  B.pseudomallei 

138  A  single  colony  of  B.  pseudomallei  grown  on  blood  agar  from  patients  sample  was  re- 

139  suspended  in  ultrapure  water.  The  suspension  was  heated  at  95  °C  for  10  min  and 

140  centrifuged  at  13500  x  g  to  pellet  the  eell  debris.  The  supernatant  was  used  as  the 

141  template  for  all  subsequent  PCR  assays.  Real  time  PCR  assay  was  done  for  gene 

142  targets  of  the  IpxO,  YLF  and  BTFC  gene  elusters  using  the  primers  and  methods 

143  deseribed  previously  (18,  19). 

144  IHA  antibody  testing 
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145  Antibody  testing  against  B. pseudomallei  antigen  was  performed  using  an  in-house 

146  method  adapted  from  Alexander  et  al,  1970  (20).  Antigen  was  prepared  from  heat 

147  killed  culture  supernatant  of  a  Sri  Lankan  B. pseudomallei,  strain  BPs7.  A  1/80  diluted 

148  antigen  preparation  was  used  to  sensitize  sheep  erythrocytes.  Serum  samples  were 

149  heat  inactivated  at  56°C  for  30  mins  and  tested  by  serial  dilution  from  1/10  to 

150  1/10,240  with  sensitized  sheep  erythrocytes  and  the  reciprocal  of  the  highest  dilution 

151  at  which  hemagglutination  occurred  was  recorded  as  end  point  titer  (20). 

152  Sample  collection  and  processing 

153  10  ml  of  whole  blood  was  collected  from  patients/volunteers  after  written  informed 

154  consent,  of  which  7ml  were  collected  into  B.D  vacutainer  mononuclear  cell 

155  preparation  tubes  (catalog  362761)  for  lymphocyte  purification.  The  lymphocytes 

156  were  purified  using  the  Ficoll  fractionation  method  as  per  manufacturer’s  instructions 

157  and  lysed  with  RLT  buffer  (Qiagen  RNeasy  mini  kit-catalog  74104),  homogenized 

158  and  stored  at  -80°C  for  total  RNA  extraction. 

159  Total  RNA  extraction  and  cDNA  synthesis 

160  Total  RNA  was  extracted  from  the  stored  cell  lysate  samples  using  the  Qiagen 

161  RNeasy  mini  kit  (catalog  no:74104)  as  per  manufacturer  protocol.  RNA  extracted 

162  from  0.6  million  PBMC’s  was  used  for  cDNA  synthesis  as  the  standard  for  all 

163  samples  analyzed  by  RT-qPCR.  cDNA  was  synthesized  using  Qiagen  First  strand  kit 

164  (catlog330401)  as  per  the  manufacturer  recommended  instructions.  The  synthesized 

165  cDNA  samples  were  stored  at  -20°C  until  further  use. 

166  Real  Time  qPCR  and  gene  expression  analysis 
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'y 

167  Qiagen  human  common  cytokines  RT  Profiler  PCR  array  (catlog  PAHS-021Z)  was 

168  used  for  this  study.  The  PCR  reaction  and  thermal  profile  recommended  by  the 

169  manufacturer  were  followed.  26  melioidosis  cases  (identified  as  confirmed  or 

170  probable  cases),  10  other  bacterial  sepsis  cases  (negative  for  B. pseudomallei)  and  5 

171  healthy  negative  controls  were  analyzed  by  RT-qPCR. 

172  Data  Analysis 

173  The  relative  gene  expression  ratio,  for  measuring  the  change  in  expression  level  of  a 

174  gene  was  calculated  by  delta  delta  CT  method  (21)  as  per  manufacturer 

175  recommendations.  The  data  was  normalized  using  actin  beta  as  the  reference 

176  housekeeping  gene.  Statistical  analysis  was  done  by  Welch’s  T-test  using  SAS  PROC 

177  MIXED,  version  9.4.  P<0.05  was  considered  as  statistically  significant. 

178 

179  RESULTS  AND  DISCUSSION 

180  A  total  of  26  cases  of  melioidosis  were  analyzed  of  which  23  were  confirmed  cases 

181  (culture  positive),  and  3  were  probable  cases  (high  antibody  titre  positive).  A  majority 

182  (n=23)  of  melioidosis  cases  had  associated  comorbidities,  and  diabetes  was  the  most 

183  common  comorbidity  (n=17)  in  this  study.  Out  of  23  confirmed  cases  of  melioidosis 

184  16  were  classified  as  septicaemic  or  bacteriaemic. 

185  The  differential  expression  pattern  of  interleukins  (IL),  interferons  (IFN),  tumor 

186  necrosis  factor  (TNF)  super  family,  transforming  growth  factor  (TGF)  super  family 

187  and  other  growth  factors  was  significant  in  melioidosis  patients  compared  to  other 

188  bacterial  sepsis  infection  cases  and  healthy  controls.  Adiponectin,  CIQ  and  collagen 

189  domain  containing  (ADIPOQ)  and  family  with  sequence  similarity  3,  member  B 
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190  (FAM3B)  were  significantly  down  regulated  in  other  bacterial  sepsis  infection  cases 

191  compared  to  healthy  controls  (fig  1). 

192  Gene  expression  profile  of  interleukins 

193  Our  study  reveals  up  regulated  expression  of  ILIO,  ILIB,  ILIRN  (interleukin  1 

194  receptor  antagonist),  IL27  and  IL8  in  melioidosis  patients  compared  to  healthy 

195  controls  (Table  1,  Fig  l).This  is  in  agreement  with  a  study  by  Weirsinga  et  al,  2007 

196  reporting  increased  mRNA  expression  of  inflammatory  response  genes  such  as  ILip, 

197  IL6,  IL15,  ILIO,  IL4,  IFNy  and  TNFa  in  melioidosis  patients  when  compared  to 

198  healthy  controls  (7).  In  our  study,  ILI6,  ILI7A,  IL23A  and  IL24  were  down  regulated 

199  while  ILIO,  ILIB  and  IL8  were  up  regulated  in  other  bacterial  sepsis  infection  cases 

200  compared  to  healthy  controls  (Fig  1).  IL16,  IL17A,  IL17B,  ILIRN,  IL22,  IL23A, 

201  IL24,  IL27,  IL3  and  IL4  were  all  up  regulated  in  melioidosis  patients  compared  to 

202  other  bacterial  sepsis  infection  cases  (Table  2,  Fig  1).  Particularly,  IL17A,  IL3  and 

203  IL4  showed  high  levels  of  gene  expression.  Previously,  expression  profiling  of 

204  interleukins  in  response  to  B. pseudomallei  infection  has  been  extensively  studied  in 

205  several  animal  models,  showing  upregulated  expression  of  interleukins  such  as  ILip, 

206  IL6,  ILIO  and  ILI2  within  72  hours  of  infection  (9-13,  15).  Elevated  levels  of 

207  expression  of  IL6,  IL8,  IL12,  IL15  and  IL18  was  also  observed  in  the  plasma  of 

208  melioidosis  patients  (15,  16). 

209  IL17A,  a  pro-inflammatory  cytokine  which  mediates  inflammatory  responses  and 

210  induces  production  of  other  cytokines,  is  particularly  expressed  at  very  high  levels  in 

211  melioidosis  patients  (including  septicaemic  and  diabetic  cohorts)  compared  to  other 

212  sepsis  infections  (Fig  1-  2).  Additionally,  IL-22  a  widely  regarded  Thl7  cytokine,  also 

213  shows  upregulated  expression  in  melioidosis  patients  compared  to  other  sepsis  cases. 
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214  IL17  and  other  Thl7  cytokines  are  linked  to  the  response  against  extra-cellular 

215  bacteria,  pathogenesis  of  diverse  autoimmune  and  inflammatory  diseases,  as  their 

216  dysregulated  expression  can  lead  to  uncontrolled  inflammatory  responses  (22,  23).  IL- 

217  17  is  also  implicated  in  excessive  tissue  damage  by  stimulating  the  production  of 

218  many  other  cytokines  including  granulocyte-colony  stimulating  factor  (G-CSF), 

219  granulocyte-macrophage  colony  stimulating  factor  (GM-CSF),  TGF-P,  TNF-a,  thus 

220  contributing  to  inflammatory  pathology  (8).  IL23,  a  key  mediator  of  inflammation  has 

221  also  been  reported  to  show  upregulated  mRNA  expression  during  B. pseudomallei 

222  infection,  implicating  its  role  in  pathogenic  host  immune  responses  (24).  Anti-IL17 

223  and  anti-IL23  therapeutic  agents  have  shown  to  be  effective  in  several  immune- 

224  mediated  inflammatory  diseases  (23,  25).  IL27,  implicated  in  regulating  B  and  T  cell 

225  activity,  has  been  reported  to  be  significantly  elevated  in  melioidosis  patients 

226  compared  to  healthy  controls  and  over  production  of  IL27  plays  a  major  role  in 

227  pathogenesis  of  sepsis  and  shock  (26).  IL27  has  also  been  identified  as  a  potential 

228  sepsis  biomarker  and  a  candidate  in  successful  therapeutic  intervention  (27,  28).As 

229  our  results  show  consistent  upregulated  expression  of  IL17,  IL23  and  IL27,  their  role 

230  in  melioidosis  disease  progression  and  therapeutic  use  should  be  further  investigated. 

231  Our  study  revealed  greater  than  3  fold  upregulation  of  IL4,  IL13,  IL17A,  IL17B, 

232  IL22,  IL23A,  IL24  and  IL27  in  the  diabetic  melioidosis  cohort  (n=17)  compared  to 

233  other  bacterial  sepsis  cases  (Fig  2). Diabetes,  a  risk  factor  for  infectious  diseases,  may 

234  play  a  role  in  neutrophil  and  T-cell  dysfunction,  possibly  mediated  by  altered  glucose 

235  metabolism  and  oxidative  stress  (29).  Studies  on  diabetic  cohorts  (mice  and  human  ) 

236  of  melioidosis  infection  shows  excessive  neutrophil  infiltration  and  impaired 

237  inflammatory  and  Thl  cytokine  responses,  leading  to  increased  susceptibility  of 

238  diabetic  individuals  to  melioidosis  (10,  30).  IL4,  a  key  regulator  of  humoral  and 


TR-17-135 

Distribution  Statement  A:  Approved  for  pubiic  reiease;  distribution  is  uniimited. 


239  adaptive  immunity,  functions  as  an  anti-inflammatory  cytokine  decreasing  production 

240  of  Thl  cells  and  related  pro-inflammatory  cytokines.  Our  findings  show  upregulation 

241  of  IL4  and  closely  related  anti-inflammatory  cytokine  IL13,  in  the  melioidosis  cohort 

242  compared  to  other  sepsis  cases,  which  is  suggestive  of  inflammatory  responses  being 

243  dysregulated.  Upregulated  IL4  expression  has  been  reported  in  melioidosis  patients 

244  and  acute  melioidosis  animal  models  (7,  9). 

245  Our  findings  also  show  significant  upregulation  of  IL17A,  IL17B,  IL23A  whereas 

246  ILIA,  ILIB  and  IL8  were  down  regulated  in  septicaemic  melioidosis  cohort  (n=16) 

247  compared  to  other  sepsis  cases  (Fig  1).  Downregulation  of  ILIA,  ILIB,  IL6,  IL8  and 

248  IL21  in  the  early  acute  phase  melioidosis  cohort  (<15  days  fever/clinical  symptoms) 

249  compared  to  other  sepsis  cases  (Fig  3)  was  also  seen,  indicating  ILIA,  ILIB  and  IL8 

250  as  potential  markers  during  the  early  stages  of  inflammation  and  being  correlated  with 

251  disease  severity.  A  study  using  a  human  lung  epithelial  cell  line  showed  that  IL8 

252  production  upon  B. pseudomallei  infection  was  lower  than  cells  infected  with  other 

253  gram  negative  bacteria  which  correlates  with  our  findings  (14).  Increased  level  of 

254  plasma  IL6  and  IL8  concentration,  being  associated  with  disease  severity  and 

255  mortality  have  also  been  reported  (15).  Immunosuppression  in  melioidosis  patients 

256  correlating  to  mortality,  associated  with  up  regulated  interleukin- IR-associated- 

257  kinase-M  expression,  leading  to  a  strong  decrease  in  capacity  to  release  pro- 

258  inflammatory  cytokines  such  as  ILIB,  TNFa  and  IL8,  after  ex-vivo  stimulation  with 

259  LPS  or  B. pseudomallei,  has  been  reported  (31).  Downregulation  of  ILIB  upon 

260  B. pseudomallei  infection  compared  with  avirulent  B.thailandensis  in  infected  lung 

261  epithelial  cells  has  also  been  reported,  suggesting  host  response  evasion  (32).  In  our 

262  findings,  we  also  see  an  upregulation  of  ILRN,  a  natural  inhibitor  of  pro- 

263  inflammatory  effects  of  ILIA,  ILIB,  in  melioidosis  patients  compared  to  other 
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264  bacterial  sepsis  cases.  Thus  ILl  and  IL8  whieh  are  key  mediators  of  inflammation  and 

265  early  innate  immune  responses,  may  serve  as  eandidate  early  diagnostic  markers  and 

266  indicators  of  disease  severity. 

267  Our  findings  reveal  an  upregulated  Th2,  Thl7  eytokine  response  and  a  down 

268  regulated  Thl  cytokine  response,  with  associated  comorbidities  such  as  diabetes 

269  playing  a  key  role  in  pathogenesis  and  severity  through  dysregulated  eytokine 

270  responses. 

271  Gene  expression  profile  of  interferons 

272  Interferon  A5  (IFNA5)  was  down  regulated  in  melioidosis  patients  compared  to 

273  healthy  eontrols  (Table  1,  Fig  1).  Interferon  BI  (IFNBI)  was  significantly  down 

274  regulated  in  other  sepsis  infeetion  eases  eompared  to  healthy  eontrols  (Fig  1). 

275  Interferon  AI(IFNAI)  and  IFNBI  showed  upregulation  whereas  IFNA5  was  down 

276  regulated  in  the  melioidosis  patients  compared  to  other  sepsis  infeetion  cases  (Table 

277  2,  Fig  1). 

278  Elevated  expression  of  interferon  gamma  (IFNy),  a  pro-inflammatory  cytokine,  has 

279  been  reported  in  human  host  and  animal  models  of  B.  pseudomallei  infection  during 

280  the  early  stages  (7,  11,  12,  16).  Our  findings  did  not  show  any  significant  upregulation 

281  of  INFy  in  melioidosis  cases  compared  to  healthy  controls  or  other  baeterial  sepsis 

282  cases,  possibly  due  to  our  samples  being  collected  at  latter  stages  while  undergoing 

283  antibiotie  treatment.  However  we  did  see  a  signifieant  upregulation  in  the  diabetic 

284  melioidosis  cohort  compared  to  other  sepsis  cases  of  infection  (Fig  2).  Interferon 

285  mediated  responses  have  been  reported  as  the  most  dominant  pathway,  with  class  I 

286  and  II  interferons  being  prominent  in  melioidosis  and  tuberculosis  infections  (33). 

287  Our  study  shows  an  upregulated  expression  of  INFAl  and  INFBl  in  melioidosis 
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288  patients  (including  diabetic  and  septicaemic  cohorts)  compared  to  other  sepsis 

289  infections  (Fig  1,2).  Interferons  a  and  P,  both  belonging  to  class  I  interferons,  play  a 

290  major  role  in  innate  immune  responses.  Dysregulated  type  I IFN  production  results  in 

291  a  damaging  cascade  of  cell  death,  inflammation,  and  immunological  host  responses 

292  that  can  lead  to  tissue  injury  and  disease  progression  (34).  Studies  have  shown  type  I 

293  IFN  responses  as  a  striking  characteristic  of  TB  infection  and  that  lack  of 

294  development  of  Thl  immunity  in  response  to  M.  tuberculosis  appears  to  be  associated 

295  with  increased  induction  of  type  1  IFNs,  leading  to  better  bacterial  survival  and  host 

296  evasion  (35).  Furthermore  a  study  also  reported  that  type  1  IFNs  suppresses  IL-I 

297  production,  providing  cellular  basis  for  the  anti-inflammatory  effects,  as  well  as  pro- 

298  bacterial  functions  of  type  I  IFNs  during  M.  tuberculosis  infection  (36).  Our  findings 

299  show  a  similar  response,  as  we  see  a  dominant  type  I  IFN  production  and  a  fairly 

300  submissive  IFNy  and  related  Thl  cytokines  production  in  melioidosis  patients.  This 

301  data  supports  our  findings  of  significant  downregulated  expression  of  ILIA,  ILIB  in 

302  melioidosis  patients  compared  to  other  sepsis  cases.  Type  I  interferons  are  also 

303  considered  as  mediators  of  endotoxic  shock  and  sepsis  induced  by  gram  negative 

304  bacteria,  with  IFNB  and  IFNARl  deliberated  as  therapeutic  targets  (37,  38).  Thus 

305  further  investigation  is  required  to  understand  the  expression  of  class  I  interferons  in 

306  relation  to  pathogenesis  of  melioidosis  and  its  role  in  diagnostic  and  therapeutic 

307  intervention. 

308  Gene  expression  profile  of  TNF  superfamily 

309  Our  findings  reveal  an  upregulation  of  TNFa  an  important  pro-inflammatory 

310  cytokine,  in  melioidosis  cases  compared  to  healthy  controls  (Table  1,  Fig  1).  Several 

311  studies  have  reported  the  upregulated  expression  of  TNFa  during  melioidosis  in 
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312  human  host  and  animal  models  of  infection  (7,  11,  12,  39).  Elevated  plasma 

313  concentrations  of  TNFa  have  been  correlated  with  disease  severity  and  mortality  in 

314  septicaemic  melioidosis  patients  (40). 

315  CD40  ligand  (CD40LG);  which  plays  a  major  role  in  B-cell  activation  and 

316  development  and  pro-inflammatory  cytokines  and  Lymphotoxin  alpha  (LTA)  were 

317  down  regulated  in  the  melioidosis  cohort  compared  to  healthy  controls  (Table  1,  fig 

318  1).  CD40L  has  been  considered  as  an  important  mediator  of  sepsis,  implicated  in 

319  platelet-mediated  activation  and  accumulation  of  neutrophils  during  inflammation 

320  (41, 42). 

321  Tumor  necrosis  factor  (ligand)  superfamily  14  (TNFSF14);  which  plays  a  major  role  in 

322  T-cell  proliferation.  Tumor  necrosis  factor  (ligand)  superfamily  4  (TNFSF4);  which  is 

323  responsible  for  Th2  cell  differentiation,  and  Tumor  necrosis  factor  (ligand)  superfamily  8 

324  (TNFSF8);  implicated  in  blocking  Thl  responses  were  up  regulated  in  the  melioidosis 

325  cohort  compared  to  other  sepsis  infection  cases  (Table  1,  Fig  1).  TNFSF4  was 

326  consistently  upregulated  in  melioidosis  patients  compared  to  other  sepsis  cases 

327  irrespective  of  factors  like  duration  of  clinical  symptoms,  antibiotic  treatment  and 

328  comorbidities  such  as  diabetes  (fig  1-3).  Upregulated  expression  of  TNFSF4  (also 

329  known  as  OX40F),  has  been  observed  in  cases  of  polymicrobial  sepsis  and 

330  autoimmune  disease  and  has  been  correlated  to  disease  severity  and  mortality  (43). 

331  Studies  have  also  shown  that  upregulated  expression  of  TNFSF4  promoted  T  cell 

332  proliferation,  increased  expression  in  CD4-I-  T  cells  and  production  of  Th2  cytokines 

333  such  as  IF4  (44,  45).  It  has  also  been  postulated  as  a  specific  biomarker  in  therapeutic 

334  interventions  for  treatment  of  sepsis/  septic  shock  and  other  autoimmune  diseases 

335  (43).  Tumor  necrosis  factor  (ligand)  superfamilyl4  (TNFSF14),  otherwise  known  as 

336  FIGHT  plays  a  major  role  in  systemic  immune  response,  particularly  in  long  term 
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337  survival  of  memory  Thl  and  Th2  cells  (46).  Tumor  necrosis  factor  (ligand)  superfamily  8 

338  (TNFSF8),  or  CD30L,  is  reportedly  expressed  in  Th2  cells  and  suppresses  Thl 

339  responses  (46).  These  findings  once  again  suggest  an  inclination  for  dominant  Th2 

340  responses  during  the  disease  progression  of  melioidosis. 

341  Gene  expression  profile  of  TGFp  superfamily 

342  Bone  morphogenetic  protein  6  (BMP6),  inhibin  beta  A  (INHBA)  and  transforming 

343  growth  factor  beta  l(TGFBl)  showed  significant  upregulation  in  melioidosis  patients 

344  compared  to  healthy  individuals  (Table  1,  Fig  1).  BMPS  was  down  regulated  in  other 

345  sepsis  infection  cases  compared  to  healthy  controls  (Fig  1).  BMPS,  BMP4,  BMP6, 

346  growth  differentiation  factor  2  (GDF2),  INHBA,  TGFBl  expressed  up  regulated 

347  expression  in  melioidosis  patients  compared  to  other  sepsis  cases,  with  BMP4  and 

348  GDF2  showing  high  level  of  gene  expression  (Table  2  Fig  1). 

349  High  level  expression  of  BMPS,  BMP6,  TGFBl  and  TGFB2  was  observed  in 

350  septicaemic  melioidosis  cohort  compared  to  other  bacterial  sepsis  cases  (Fig  1). 

351  TGFB2  was  expressed  at  high  level  in  the  early  acute  phase  (<15  days  of 

352  fever/clinical  symptoms)  melioidosis  cohort  compared  to  sepsis  control  (Fig  3). 

353  TGFp  was  upregulated  during  melioidosis  infection,  with  increased  levels  being 

354  correlated  to  severe  melioidosis  in  human  hosts  (47).  Our  study  revealed  a  consistent 

355  upregulation  of  TGFBl  in  melioidosis  patients  compared  to  other  sepsis  cases 

356  irrespective  of  factors  like  duration  of  clinical  symptoms,  antibiotic  treatment  and 

357  comorbidities  such  as  diabetes  (fig  1-3).  An  experimental  murine  model  of 

358  melioidosis,  revealed  that  an  inhibition  of  TGF-P  with  a  selective  TGF-P  antibody 

359  had  a  protective  effect,  with  reduction  in  inflammation,  reduced  bacterial  load  and 

360  organ  damage,  thus  indicating  the  role  of  TGF-P  in  pathogenesis  of  melioidosis  (47). 
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361  Several  other  studies  have  also  shown  the  crucial  role  of  TGF-P  in  immune 

362  regulation,  where  it  induces  Foxp3,  a  master  regulator  of  Tregs  in  naive  T  cells,  with 

363  suppression  of  pro-inflammatory  cytokines  such  as  IFNy  and  enhanced  production  of 

364  anti-inflammatory  cytokines  (48,  49).  It  has  also  been  identified  as  an  inducer  of  T17 

365  cell  differentiation  (48,  50).  These  studies  further  support  our  findings  of  increased 

366  Thl7  cytokine  production  and  suppression  of  Thl  cytokines  in  melioidosis  patients. 

367  BMPS  and  BMP6  were  consistently  up  regulated  in  melioidosis  patients  compared  to 

368  other  sepsis  cases  irrespective  of  factors  like  duration  of  clinical  symptoms,  antibiotic 

369  treatment  and  comorbidities  such  as  diabetes  (Fig  1-3).  BMP’s  play  a  major  role  in 

370  formation  and  repair  of  bone  and  cartilage,  cell  proliferation,  differentiation  and 

371  apoptosis  (51).  INHBA  over  expression  has  been  associated  with  increased  cell 

372  proliferation  and  poor  disease  outcome  in  several  types  of  carcinomas  (52,  53). 

373  Further  studies  are  needed  to  elucidate  the  mechanisms  of  BMP  signaling  pathways 

374  and  INHBA  expression  in  relation  to  pathogenesis  of  melioidosis. 

375  Gene  expression  profile  of  Growth  factors 

376  Platelet-derived  growth  factor  alpha  polypeptide  (PDGFA)was  upregulated 

377  significantly  in  melioidosis  patients  compared  to  healthy  individuals  (Table  1,  Fig  1). 

378  PDGFA,  thrombopoietin  (THPO),  ciliary  neurotrophic  factor  (CNTF),  macrophage 

379  colony  stimulating  factor  (M-CSF  or  CSFl),  C-fos  induced  growth  factor  or  vascular 

380  endothelial  growth  factor  D  (FIGF)  showed  upregulation  in  melioidosis  cases 

381  compared  to  other  sepsis  infection  cases  (Table  2,  Fig  1).  A  down  regulated 

382  expression  in  THPO  was  observed  in  other  sepsis  cases  compared  to  healthy 

383  individuals  (Fig  1). 
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384  PDGF  is  an  important  growth  factor  that  plays  a  crucial  role  in  blood  vessel  formation 

385  (angiogenesis)  and  regulates  cell  growth  and  differentiation.  THPO  stimulates  the 

386  production  and  differentiation  of  megakaryocytes,  thus  regulating  platelet  production. 

387  FIGF  plays  an  active  role  in  angiogenesis  and  vascular  endothelial  cell  growth  (54). 

388  Increased  expression  of  PDGF  is  seen  in  severe  bacterial  infections,  implicating  the 

389  role  of  angiogenic  factors  in  endothelial  dysfunction  leading  to  disease  pathogenesis 

390  (54).  PDGF  has  also  been  suggested  as  a  biomarker  of  sepsis,  related  to  vascular 

391  endothelial  damage  (55).  Our  findings  also  agree  with  these  reports  as  we  see  an 

392  increased  expression  of  growth  factors,  which  play  a  role  in  endothelial  function. 

393  A  down  regulated  expression  of  leukemia  inhibitory  factor  (LIF);  an  IL6  class 

394  cytokine  that  inhibits  cell  differentiation  and  a  similar  cytokine  OSM  (oncostatin  M), 

395  was  observed  in  the  early  acute  phase  melioidosis  cases  (<15  days  fever/clinical 

396  symptoms)  compared  to  other  sepsis  cases  (Fig  3).  M-CSF  or  CSFl,  FIGF  and 

397  PDGFA  were  consistently  up-regulated  in  melioidosis  patients  compared  to  other 

398  sepsis  cases  irrespective  of  factors  like  duration  of  clinical  symptoms,  antibiotic 

399  treatment  and  comorbidities  such  as  diabetes  (Fig  1-3).  Studies  with  experimental 

400  mice  models  of  melioidosis  have  revealed  a  upregulation  of  mRNA  for  macrophage 

401  colony  stimulating  factor  (CSFl  or  M-CSF),  granulocyte  macrophage  colony 

402  stimulating  factor  (CSF2  or  GM-CSF),  granulocyte  colony  stimulating  factor  (CSF3 

403  or  G-CSF)  at  day  3  post  infection,  correlating  with  peak  bacterial  load  and  extensive 

404  infiltration  of  leucocytes  (56).  Colony  stimulating  factors  are  glycoproteins,  necessary 

405  for  the  survival,  proliferation  and  differentiation  of  hematopoietic  progenitor  cells  of 

406  the  myeloid  and  erythroid  lineage.  M-CSF  enhances  the  survival  and  activation  of 

407  cells  of  the  monocyte  lineage,  while  GM-CSF  and  G-CSF  increases  accumulation 

408  and  activation  of  both  neutrophils  and  macrophages  (56). While  colony- stimulating 
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409  factors  play  a  crucial  role  in  innate  immune  responses  and  host  defense,  their  high 

410  level  of  expression  during  melioidosis  may  instead  contribute  to  disease  pathogenesis. 

411  Limitations  of  the  study 

412  The  main  limitation  of  our  study  was  that  the  melioidosis  patient  samples  were 

413  collected  after  start  of  antibiotic  treatment  which  may  affect  immunocompetant  cells, 

414  which  in  turn  affects  the  cytokine  profiles  studied  here.  Studies  have  shown  that 

415  antibiotics  like  meropeneum  exert  an  immunomodulatory  effect,  affecting  the 

416  production  of  some  cytokines  in  PBMC’s  (57).  This  may  have  been  the  main  reason, 

417  as  to  why  we  could  not  see  any  significant  differential  expression  of  some  key 

418  inflammatory  response  cytokines  such  as  IFNy.  Duration  of  clinical  symptoms  ranged 

419  from  >10  days  to  >90  days  and  duration  of  antibiotics  treatment  ranged  from  3  days 

420  to  >30  days  at  the  time  of  blood  collection  for  all  the  melioidosis  samples.  Since  our 

421  sample  collection  was  nationwide,  duration  between  patient  identification/disease 

422  confirmation  and  sampling  was  substantial  due  to  logistical  issues.  Thus,  due  to 

423  varying  and  the  wide  range  of  acute  phase  in  each  of  the  samples  analyzed  and  less 

424  number  of  samples  with  <15  days  of  fever/clinical  symptoms  duration  (n=5)  we  could 

425  not  see  any  statistically  significant  differential  expression  of  some  of  the 

426  inflammatory  response  genes  involved  in  early  innate  immune  responses.  However, 

427  our  results  showed  a  consistent  up  regulated  expression  of  interleukins;  IL4,  IL17A, 

428  IL23A,  IL24,  interferons;  IFNAl,  IFNBl,  TNF  superfamily;  TNFSF4  (OX40L),  TGF 

429  superfamily;  BMP3,  BMP6,  TGFBl;  other  growth  factors;  CSFl,  FIGF  and  PDGFA 

430  in  melioidosis  patients  compared  to  other  sepsis  cases,  irrespective  of  comorbidities, 

431  duration  of  fever/clinical  symptoms  and  antibiotic  treatment,  indicating  their 
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432  differential  expression  during  melioidosis  infection. Our  findings  suggest  a  domination 

433  of  Th2  and  Thl7  type  responses  during  disease  pathogenesis  of  melioidosis. 

434  As  diabetes  was  seen  as  a  major  comorbidity  in  our  experimental  cohort,  we  analyzed 

435  our  data  to  see  if  there  was  any  significant  differential  expression  between  diabetic 

436  melioidosis  cases  and  non-diabetic  melioidosis  cases.  The  gene  expression  pattern 

437  between  these  two  groups  were  comparable  and  we  could  not  find  any  statistically 

438  significant  differential  expression,  indicating  that  the  differential  expression  was 

439  largely  due  to  melioidosis  infection  (Tables  SI  and  S2). 

440  Conclusion 

441  Our  study  revealed  differential  gene  expression  of  key  cytokines  involved  in  human 

442  host  responses  that  can  distinguish  melioidosis  cases  from  sepsis  infections  caused  by 

443  other  pathogens  and  healthy  individuals.  Low  level  of  expression  of  key  inflammatory 

444  mediators;  ILIA,  IL IB  and  IL8  were  seen  in  melioidosis  patients  in  early  acute  phase 

445  and  with  septicaemia  compared  to  other  sepsis  infection  cases.  These  findings 

446  indicate  that  differentially  expressed  genes  should  be  validated  during  different  stages 

447  of  infection  for  their  potential  as  disease  biomarkers  for  diagnostic  purposes  and 

448  monitoring  disease  progression.  Our  results  also  show  an  elevated  expression  of  Thl7 

449  cytokines  such  as  IL17,  IL22  and  TGF  P  which  act  as  an  inducer  of  Thl7  cytokines. 

450  Th2  cytokines  such  as  IL3,  IL4  and  IL13  were  also  upregulated  along  with  type  I 

451  interferons  and  TNFSF  cytokines,  which  are  known  to  be  inducer’s  of  Th2  cytokines 

452  and  suppressors  of  Thl  responses.  These  results  may  indicate  a  dominant  Th2  and 

453  Thl7  type  cytokine  responses,  suggesting  that  their  dysregulation  may  play  an 

454  important  role  in  disease  pathogenesis  and  progression.  IL17,  IL23  and  IL27,  already 

455  implicated  in  therapeutic  intervention  of  several  inflammatory  diseases  should  be 
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456  further  investigated  for  their  role  in  disease  progression  and  therapeutic  approaches  in 

457  melioidosis. 

458  Our  future  studies  shall  be  aimed  at  studying  gene  expression  profiles  in  early  and  late 

459  acute  phases  of  melioidosis  to  evaluate  candidate  genes  which  can  serve  as  disease 

460  and  diagnostic  biomarkers  in  different  stages  of  infection.  Based  on  these  biomarkers 

461  if  the  antibiotic  treatment  regime  can  be  adjusted  it  would  bring  benefits  to  the 

462  patients  by  reducing  the  hospital  stay.  We  would  expand  our  studies  further,  with  a 

463  larger  sample  size  in  each  category  of  sample  type,  focusing  on  specific  immune 

464  response  genes  showing  differential  expression,  to  further  understand  their  role  in 

465  disease  pathogenesis,  susceptibility  and  severity  associated  with  major  comorbidities 

466  such  as  diabetes. 

467 
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682  TABLES 

683  Table  1.  Cytokines  showing  significant  differential  expression  in  PBMC’s  of 

684  melioidosis  patients  (n=26)  compared  to  healthy  negative  controls  (n=5) 
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701  confidence  interval;  Relative  expression  ratio  >1.5  indicates  upregulation  and  <0.5 

702  indicates  downregulation  in  the  experimental  group  compared  to  control  group;  Gene 

703  targets  showing  significant  differential  expression  where  P-value  <0.05 


Gene 

Gene  description 

Relative  expression  ratio 
[95%  Cl] 

p  bbb 

value 

ILIB 

Interleukin  1  beta 

2.504  [1.229  -5.100] 

0.(^ 

ILIRN 

Interleukin  1  receptor 
antagonist 

1.62  [1.023  -2.564] 

0.(M)^ 

688 

IL8 

Interleukin  8 

2.953  [1.394-6.257] 

0.0062 

689 

ILIO 

Interleukin  10 

2.257  [1.180-4.319] 

0.0158 

600 

IL27 

Interleukin  27 

4.022  [1.632-9.915] 

0.0039 

INFA5 

Interferon  alpha  5 

0.189  [0.037-0.961] 

0.0A34 

TNF 

Tumour  necrosis  factor 

2.248  [1.082-4.670] 

CD40LG 

CD40  ligand 

0.502  [0.336-0.752] 

LTA 

Lymphotoxin  alpha 

0.327  [0.190-0.565] 

O.QEBOS 

BMP6 

Bone  morphogenetic 
protein  6 

2.946  [1.470-5.904] 

o-Ws 
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Inhibin  beta  A 

6.07  [2.652  -  13.891] 
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TGFBl 

Transforming  growth 
factor  1 
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0.(^6 

698 

PDGFA 

Platelet  derived  growth 
factor  alpha  polypeptide 

2.86  [1.444-5.667] 

0.0066 
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705  Table  2.  Cytokines  showing  significant  differential  expression  in  PBMC’s  of 

706  melioidosis  patients  (n=26)  compared  to  other  sepsis  cases  (n=10) 


Gene 

symbol 

Gene  description 

Relative  expression  ratio 
[95%  Cl] 

P  value 

IL3 

Interleukin  3 

9.38  [1.773  -49.626] 

0.0107 

IL4 

Interleukin  4 

6.024  [1.153- 31.479] 

0.0344 

IL16 

Interleukin  16 

1.896  [1.152-3.121] 

0.0157 

IL17A 

Interleukin  17  alpha 

16.32  [3.193  -  83.421] 

0.0017 

IL17B 

Interleukin  17  beta 

2.939  [1.486-5.811] 

0.003 

ILIRN 

Interleukin  1  receptor 
antagonist 

1.747  [1.091-2.796] 

0.0216 

IL22 

Interleukin  22 

3.022  [1.207-7.565] 

0.0206 

IL23A 

Interleukin  23  alpha 

2.792  [1.329-5.866] 

0.0092 

IL24 

Interleukin  24 

2.991  [1.240  -7.214] 

0.0173 

IL27 

Interleukin  27 

3.089  [1.203  -7.932] 

0.0206 

INFAl 

Interferon  alpha  1 

4.034  [1.358  -11.984] 

0.014 

INFA5 

Interferon  alpha  5 

0.2  [0.057  -  0.704] 

0.0152 

INFBl 

Interferon  beta  1 

3.206  [1.056-9.735] 

0.0407 

TNFSF4 

Tumour  necrosis  factor  super 
family  4 

2.349  [1.167  -4.728] 

0.0202 

TNFSF8 

Tumour  necrosis  factor  super 
family  8 

1.606  [1.004  -2.571] 

0.0484 

TNFSF14 

Tumour  necrosis  factor  super 
family  14 

2.353  [1.171  -4.728] 

0.0186 

BMP3 

Bone  morphogenetic  protein  3 

5.305  [2.319  -  12.135] 

0.0003 

BMP4 

Bone  morphogenetic  protein  4 

18.765  [1.479  -238.054] 

0.0271 

BMP6 

Bone  morphogenetic  protein  6 

2.776  [1.214 -6.344] 

0.0192 

GDF2 

Growth  differentiation  factor  2 

11.112  [1.105  -111.704] 

0.0421 

INHBA 

Inhibin  beta  A 

4.635  [1.205  -17.822] 

0.0282 
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707 

708 

709 

710  PB 

711  MC, 

712  perip 

713  heral 

714  blood  mononuclear  cells;  Cl,  confidence  interval;  Relative  expression  ratio  >1.5 

715  indicates  upregulation  and  <0.5  indicates  downregulation  in  the  experimental  group 

716  compared  to  control  group;  Gene  targets  showing  significant  differential  expression 

717  where  P- v  alue  <0.05 

718 


TGFBl 

Transforming  growth  factor  1 

2.006  [1.374 -2.931] 

0.0007 

PDGFA 

Platelet  derived  growth  factor 
alpha  polypeptide 

2.317  [1.065  -5.038] 

0.0357 

THPO 

Thrombopoietin 

4.213  [1.042-  17.040] 

0.0441 

CNTF 

Ciliary  neurotrophic  factor 

2.222  [1.023  -4.829] 

0.0441 

CSFl 

Colony  stimulating  factor  1 

2.456  [1.451  -4.156] 

0.0017 

FIGF 

C-fos  induced  growth  factor 

3.912  [1.561  -9.802] 

0.0049 
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719  FIGURES 

720  Fig  1:  Relative  differential  gene  expression  of  cytokines  in  melioidosis  patients 

721  compared  to  patients  with  sepsis  infection  due  to  other  pathogens  and  healthy 

722  negative  controls.  Significant  relative  gene  expression  changes  in  PBMC’s  from 

723  melioidosis  patients  (n=26)  and  septicaemic  melioidosis  patients  (n=16)  compared  to 

724  sepsis  controls  (n=10)  and  healthy  controls  (n=5).  Expression  levels  were  normalized 

725  against  beta  actin  as  the  reference  house  keeping  gene.  Relative  expression  ratio  >1.5 

726  considered  as  upregulation  and  <0.5  was  considered  as  downregulation.  *  indicates 

727  relative  expression  ratio  is  significantly  different  (P<0.05) 

728  Fig  2:  Relative  differential  gene  expression  of  cytokines  in  diabetic  melioidosis 

729  patients  compared  to  patient  with  sepsis  infection  due  to  other  pathogens  and 

730  healthy  negative  controls.  Relative  gene  expression  in  PBMC’s  from  diabetic 

731  melioidosis  patients  (n=17),  compared  to  sepsis  controls  (n=10)  and  healthy  controls 

732  (n=5).  Expression  levels  were  normalized  against  beta  actin  as  the  reference  house 

733  keeping  gene.  Relative  expression  ratio  >1.5  was  considered  as  upregulation  and  <0.5 

734  considered  as  downregulation.  *  indicates  relative  expression  ratio  is  significantly 

735  different  (P<0.05) 

736  Fig  3:  Relative  differential  gene  expression  of  cytokines  in  melioidosis  patients 

737  with  respect  to  duration  of  fever/clinical  symptoms  and  antibiotics  treatment, 

738  compared  to  patients  with  sepsis  infection  due  to  other  pathogens.  Relative  gene 

739  expression  in  PBMC’s  from  melioidosis  patients  with  <15  days  of  fever  (n=5), 

740  melioidosis  patients  with  >15  days  of  fever  (n=21),  melioidosis  patients  with  <15 

741  days  of  treatment  with  antibiotics  (n=14),  melioidosis  patients  with  >15  days  of 

742  treatment  with  antibiotics  (n=8),  compared  to  sepsis  controls  (n=10).  Expression 
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743  levels  were  normalized  against  beta  actin  as  the  reference  house  keeping  gene. 

744  Relative  expression  ratio  >1.5  was  considered  as  upregulation  and  <0.5  considered  as 

745  downregulation.  *  indicates  relative  expression  ratio  is  significantly  different  (P<0.05) 

746 
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747  Table  1:  Primer  details*  for  Gene  expression  analysis 


Primer 

Name 

Forward  Sequence 

Reverse  Sequence 

Product 

Size 

GAPDH2 

TGACAACTTTGGTATYCGTGGAAGG 

AGGCAGGGATGATGTTCTGGAGAG 

134  bp 

18srRNA 

GCTTAATTTGACTCAACACGGGA 

AGCTATCAATCTGTCAATCCTGTC 

69  bp 

PLCEl 

GCCCAAAGCAAGTGGAAAGG 

TCTTCACCTGGGTTAAACATGC 

700  bp 

118 

CAGAGACAGCAGAGCACACA 

GGCAAAACTGCACCTTCACA 

158  bp 

MICE 

CACCCAGGCTGCAGTTCACT 

CGGGAGTCTGAGGTACGAGAA 

88  bp 

PSMB8 

GATCTCCAGAGCTCGCTTTA 

GTTCACCCGTAAGGCACTAA 

200  bp 

CCL5 

CCATGAAGGTCTCCGCGGCAC 

CCTAGCTCATCTCCAAAGAG 

361  bp 

IL6 

GGTACATCCTCGACGGCATCT 

GTGCCTCTTTGCTGCTTTCAC 

88  bp 

IL18 

GCTTGAATCTAAATTATCAGTC 

CAAATTGCATCTTATTATCATG 

335  bp 

DNMTIA 

ACCGCTTCTACTTCCTCGAGGCCTA 

GTTGCAGTCCTCTGTGAACACTGTG 

G 

335  bp 

DNMT3B 

AATGTGAATCCAGCCAGGAAAGGC 

ACTGGATTACACTCCAGGAACCGT 

191  bp 

HDACl 

CCAAGTACCACAGCGATGAC 

TGGACAGTCCTCACCAACG 

no  bp 

HDAC2 

TGAAGGAGAAGGAGGTCGAA 

GGATTTATCTTCTTCCTTAACGTCTG 

124  bp 

DNMT3A 

CGTTGGCATCCACTGTGAATGA 

TTACACACACGCAAAATACTCCTT 

551bp 

IFNy 

CCAACGCAAAGCAATACATGA 

CCTTTTTCGCTTCCCTGTTTTA 

79  bp 

TNFa 

GGA  GAA  GGG  TGA  CCG  ACT  CA 

CTG  CCC  AGA  CTC  GGC  AA 

70  bp 

ILip 

GCAAGGGCTTCAGGCAGGCCGCG 

GGTCATTCTCCTGGAAGGTCTGTGG 

GC 

96  bp 

IL4 

ACTTTGAACAGCCTCACAGAG 

TTGGAGGCAGCAAAGATGTC 

74  bp 

IL15 

GTCTTCATTTTGGGCTGTTTCAGT 

CCTCACATTCTTTGCATCCAGATTCT 

316  bp 

TLR2 

GGGTCATCATCAGCCTCTCC 

AGGTCACTGTTGCTAATGTAGGTG 

181  bp 

TLR4 

CAGAGTTGCTTTCAATGGCATC 

AGACTGTAATCAAGAACCTGGAGG 

282  bp 
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PSME2 

GGGAATGAGAAAGTCCTGTCC 

TCAATCTTGGGGATCAGGTG 

113  bp 

IL12 

CCAAGAACTTGCAGCTGAAG 

TGGGTCTATTCCGTTGTGTC 

355bp 

HDAC4 

GAGAGACTCACCCTTCCCG 

CCGGTCTGCACCAACCAAG 

240bp 

PSMA5 

AAGCCCATGTTGCTTTTGGG 

GGCGAACGGACATAGGCTAA 

112bp 

PSMB2 

AGAGGGCAGTGGAACTCCTT 

AGGTTGGCAGATTCAGGATG 

72bp 

HMGBl 

ACATCCAAAATCTTGATCAGTTA 

AGGACAGACTTTCAAAATGTTT 

122bp 

HLADMB 

ATGTGAAATCCTTTGGAGTCCCA 

GGCATCTTTACAGAGCAGAGCAT 

145bp 

748  *GAPDH  and  ISsrRNA  primers  are  the  house  keeping  genes.  PLCEl  primer  pair  tested  to 

749  amplify  700bp  genomic  region  of  PLCEl  is  used  as  genomic  DNA  control. 


750 
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751  Table  2.  mRNA  expression  in  PBMC’s  of  melioidosis  patients  (n=30)  compared  to  other 

752  sepsis  cases  (n=10)  and  healthy  negative  controls  (n=10) 

753 


Melioidosis  vs  Healthy  controls 

Melioidosis  vs  Sepsis  controls 

Sepsis  cases  vs  Healthy  controls 

Gene 

Relative  expression 
ratio  [95%  Cl] 

P 

value 

Relative  expression 
ratio  [95%  Cl] 

P 

value 

Relative  expression 
ratio  [95%  Cl] 

P 

value 

HMGBl 

0.83[0.444, 1.539] 

0.5324 

0.26[0.132, 0.510]* 

0.0005 

3.18[1.539,6.578]* 

0.0036 

IL6 

1.56[0.670, 3.623] 

0.2912 

0.89[0.360, 2.189] 

0.7867 

1.76[0.710,4.344] 

0.2079 

IL8 

1.65[0.727, 3.733] 

0.2187 

0.43[0.169, 1.080] 

0.0699 

3.86[1.411, 10.540]* 

0.0114 

ILlp 

1.33[0.554, 3.205] 

0.5005 

0.77[0.258, 2.320] 

0.6254 

1.72[0.516,5.738] 

0.3545 

IFNy 

1.34[0.826, 2.189] 

0.2174 

0.88[0.496, 1.572] 

0.6522 

1.52[0.800,2.899] 

0.1864 

TNFa 

1.22[0.705,2.124] 

0.4573 

0.81[0.408, 1.598] 

0.5173 

1.52[0.741,3.100] 

0.2365 

IL15 

1.31  [0.75 1,2.296] 

0.3254 

0.57[0.258, 1.279] 

0.1610 

2.29[1.025,5.103]* 

0.0443 

IL4 

1.78[0.115, 27.531] 

0.3533 

4.09[1. 178, 14.173]* 

0.0366 

0.44[0.065,2.942] 

0.2329 

TLR2 

1.16[0.666, 2.034] 

0.5844 

0.44[0.196,1.009] 

0.0522 

2.62[1. 182,5.803]* 

0.0212 

TLR4 

0.97[0.532, 1.777] 

0.9247 

0.37[0.141,0.974]* 

0.0448 

2.62[1.025,6.700]* 

0.0450 

MICB 

0.95[0.533, 1.699] 

0.8607 

0.33[0.183,0.582]* 

0.0006 

2.92[1.534,5.545]* 

0.0026 

HLADMB 

0.84[0.502, 1.396] 

0.4814 

0.69[0.396,1.204] 

0.1809 

1.21[0.682,2.156] 

0.4897 

PSMB2 

0.88[0.513, 1.519] 

0.6403 

0.25[0.121,0.508]* 

0.0008 

3.56[1.708,7.420]* 

0.0021 

PSME2 

1.24[0.763, 1.998] 

0.3757 

0.38[0.195,0.726]* 

0.0061 

3.28[1.684,6.389]* 

0.0017 

PSMB8 

0.83[0.456, 1.519] 

0.5356 

0.27[0.131,0.565]* 

0.0014 

3.06[1.443,6.480]* 

0.0060 

PSMA5 

0.75[0.421, 1.328] 

0.3086 

0.84[0.412,1.706] 

0.6102 

0.89[0.448,1.778] 

0.7301 

DNMTIA 

0.66[0.376, 1.171] 

0.1494 

0.51[0.211, 1.213] 

0.1171 

1.31[0.542,3.169] 

0.5224 

DNMT3A 

0.68[0.416, 1.097] 

0.1087 

0.68[0.356,1.300] 

0.2284 

0.99[0.533,1.849] 

0.9793 

DNMT3B 

1.07[0.636, 1.814] 

0.7770 

0.22[0.088,0.539]* 

0.0040 

4.94[1.948,12.503]* 

0.0031 

HDACl 

0.78[0.556, 1.095] 

0.1441 

0.50[0.353,0.719]* 

0.0006 

1.55[1.059,2.267]* 

0.0266 

HDAC2 

0.68[0.377, 1.212] 

0.1799 

0.28[0.126,0.642]* 

0.0048 

2.38[1.047,5.403]* 

0.0398 

HDAC4 

0.99[0.588, 1.652] 

0.9549 

0.64[0.322,1.270] 

0.1880 

1.54[0.786,3.025] 

0.1911 
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754 

755  PBMC,  peripheral  blood  mononuclear  cells;  Cl,  confidence  interval;  Relative  expression  ratio 

756  >1  indicates  up  regulation  and  <0.5  indicates  down  regulation  in  the  experimental  group 

757  compared  to  control  group;  *  indicates  statistically  significant  differential  expression  where 

758  P-value  <0.05,  P  values  calculated  by  paired  t-tests 

759 

760 

761 
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762 

763  Figure  1.  Relative  expression  of  genes  involved  in  immune  response  and 

764  epigenetic  regulation  in  melioidosis  patients  compared  to  patients  with  sepsis 

765  infections  caused  by  other  pathogens  and  healthy  controls.  Statistically  significant 

766  differential  expression  of  genesin  PBMC’s  from  melioidosis  patients  (n=30),  septicaemic 

767  melioidosis  patients  (n=18)  compared  to  sepsis  controls  (n=10).  Significant  differential 

768  expression  was  not  observed  among  melioidosis  patients  compared  to  healthy  controls 

769  (n=10)  while  significant  differential  expression  was  observed  among  sepsis  controls  and 

770  the  healthy  controls.  Expression  levels  were  normalized  against  18srRNA.  Relative 

771  expression  ratio  >1.5  considered  as  up  regulation  and  <0.5  considered  as  down  regulation, 

772  with  P<0.05  considered  statistically  significant. 

773 
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■  sepsis  vs  healthy  ■  melioidosis  vs  sepsis 

■  septic  melioidosis  vs  sepsis  ■  melioidosis  vs  healthy 
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786  Figure  2.  Relative  expression  of  genes  involved  in  immune  responses  and 

787  epigenetic  regulation  in  melioidosis  patients  compared  to  patients  with  sepsis 

788  infections  caused  by  other  pathogens,  in  relation  to  duration  of  fever/clinical 

789  symptoms  and  antibiotics  treatment.  Differential  gene  expression  in  PBMC’s  from 

790  melioidosis  patients  with  <15  days  of  fever  (n=4),  melioidosis  patients  with  >15  days  of 

791  fever  (n=25),  melioidosis  patients  with  <15  days  of  treatment  with  antihiotics  (n=15), 

792  melioidosis  patients  with  >15  days  of  treatment  with  antihiotics  (n=12)  compared  to  sepsis 

793  controls  (n=10),  did  not  change  due  to  the  duration  of  fever  or  duration  of  treatment 

794  with  antibiotics.  Expression  levels  were  normalized  against  18srRNA.  Relative 

795  expression  ratio  >1.5  considered  up  regulated  and  <0.5  considered  as  down  regulation,  with 

796  P<0.05  considered  statistically  significant. 

797 
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Figure  3.  Relative  expression  of  genes  involved  in  immune  responses  and 
epigenetic  regulation  in  melioidosis  patients  compared  to  patients  with  sepsis 
infections  caused  by  other  pathogens,  in  relation  to  associated  comorbidities. 

Differential  gene  expression  in  PBMC’s  from  melioidosis  patients  with  risk  factors 
(n=27),  alcoholic  melioidosis  patients  (n=8),  diabetic  melioidosis  patients  (n=20), 
compared  to  sepsis  controls  (n=10),  did  not  change  among  patientss  presented  with 
different  risk  factors  for  melioidosis.  Expression  levels  were  normalized  against 
18srRNA.  Relative  expression  ratio  >1.5  was  considered  as  up  regulated  and  <0.5  was 
considered  as  down  regulated,  with  P<0.05  considered  statistically  significant. 
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819  SI:  List  of  gene  targets  investigated 


Abbreviated  name 

Full  Name 

Biological  Role 

GAPDH 

Glyceraldehyde  3 -phosphate 
dehydrogenase 

Involved  in  several  metaholic  processes  such  as  glycolysis 

ISSrRNA 

18S  rihosomal  RNA 

Structural  RNA  and  basic  component  of  eukaryotic  cells 

PLCEl 

1  -Phosphatidylinositol-4,5- 
hisphosphate 

phosphodiesterase  epsilon- 1 

Involved  in  intracellular  responses  involving  cell  growth  and 
differentiation 

ILlp 

Interleukin  1  beta 

Pro-inflammatory  cytokine,  important  mediator  of  the 
inflammatory  response,  involved  in  a  variety  of  cellular 
activities,  including  cell  proliferation,  differentiation  and 
apoptosis. 

IL4 

Interleukin  4 

Anti-inflammatory  cytokine  that  induces  differentiation  of 
naive  helper  T  cells  (ThO  cells)  to  Th2  cells,  decreasing 
production  of  Thl  cells.  It  is  a  key  regulator  in  humoral  and 
adaptive  immunity  playing  a  major  role  in  stimulation  of 
activated  B-cell  and  T-cell  proliferation. 

IL6 

Interleukin  6 

Pro-inflammatory  cytokine  and  an  anti-inflammatory 
myokine.  It  is  secreted  by  T  cells  and  macrophages  to 
stimulate  immune  response 

IL8 

Interleukin  8 

Chemokineassociated  with  inflammation,  induces  chemotaxis 
in  target  cells,  mainly  involved  in  neutrophil  recruitment  and 
degranulation. 

ILIO 

Interleukin  10 

Anti-inflammatory  cytokine  involved  in  immune  regulation 
and  inflammation.  It  down  regulates  the  expression  of  Thl 
cytokines,  MHC  class  II  antigens,  and  co-stimulatory 
molecules  on  macrophages.  It  also  enhances  B  cell  survival, 
proliferation,  and  antibody  production.  IL-10  can  block  NF- 
kB  activity,  and  is  involved  in  the  regulation  of  the  JAK- 
STAT  signaling  pathway. 

IL12 

Interleukin  12 

Pro-inflammatory  cytokine  involved  in  the  differentiation  of 
naive  T  cells  into  Thl  cells.  It  is  involved  in  stimulation  and 
growth  T  cells  and  in  production  of  interferon-gamma  (IFN-y) 
and  tumor  necrosis  factor-alpha  (TNF-a)  from  T  cells  and 
natural  killer  (NK)  cells.  It  reduces  IF -4  mediated  suppression 
of  IFN-y. 

IL15 

Interleukin  15 

Pro-inflammatory  cytokine  which  regulates  T  and  natural 
killer  (NK)  cell  activation  and  proliferation. 
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IL18 

Interleukin  18 

Pro-inflammatory  cytokine  involved  in  inflammation  and 
cell-mediated  immunity  along  with  IL12 

CCL5 

Chemokine  (C-C  motif) 
ligand  5.  Also  known  as 
RANTES  (regulated  on 
activation,  normal  T  cell 
expressed  and  secreted). 

Chemokine,  which  is  chemotactic  for  T  cells,  eosnophils  and 
basophils,  plays  an  active  role  in  recruiting  leukocytes  into 
inflammatory  sites 

IFNy 

Interferon  gamma 

Pro-inflammatory  cytokine  that  is  critical  for  innate  and 
adaptive  immunity  against  viral,  some  bacterial  and  protozoal 
infections.  It  is  an  important  activator  of  macrophages  and 
inducer  of  Class  II  major  histocompatibility  complex  (MHC) 
molecule  expression. 

TNFa 

Tumor  necrosis  factor  alpha 

Pro-inflammatory  cytokine  involved  in  systemic 
inflammation  and  immune  regulation 

HMGBl 

High  mobility  group  hox  1 
protein,  also  known  as  high- 
mohility  group  protein  1 
(HMG-1) 

Cytokine  mediator  of  inflammation  secreted  by  activated 
macrophages  and  monocytes 

TLR2 

Toll-like  receptor  2 

Plays  a  fundamental  role  in  pathogen  recognition  and 
activation  of  innate  immunity.  This  gene  is  expressed 
abundantly  in  peripheral  blood  leukocytes  and  mediates  host 
response  to  gram-positive  bacteria. 

TLR4 

Toll-like  receptor  4 

Plays  a  fundamental  role  in  pathogen  recognition  and 
activating  the  innate  immune  system.  It  is  well-known  for 
recognizing  lipopolysaccharide  (EPS),  a  component  present 
in  many  gram-negative  bacteria,  thus  mediates  its  host 
responses 

MICE 

MHC  class  I  polypeptide- 
related  sequence  B 

Heavily  glycosylated  protein  which  is  a  ligand  for  the 

NKG2D  type  II  receptor.  Binding  of  the  ligand  activates  the 
cytolytic  response  of  natural  killer  (NK)  cells,  CD8  alpha  beta 

T  cells,  and  gamma  delta  T  cells  which  express  the  receptor. 

PSMB8 

Proteasome  suhunit  beta 
type-8  also  known  as  20S 
proteasome  subunit  beta-5i 

Eorms  a  pivotal  component  for  the  Ubiquitin-Proteasome 
System  (UPS)  involved  in  protein  ubiquitination  and 
subsequent  proteolysis  and  degradation  which  are  important 
mechanisms  in  the  regulation  of  the  cell  cycle,  cell  growth 
and  differentiation,  gene  transcription,  signal  transduction  and 
apoptosis.  During  the  antigen  processing  for  the  major 
histocompatibility  complex  (MHC)  class-I,  the  proteasome  is 
the  major  degradation  machinery  that  degrades  the  antigen 
and  present  the  resulting  peptides  to  cytotoxic  Tcells. 

TR-17-135 

Distribution  Statement  A:  Approved  for  pubiic  reiease;  distribution  is  uniimited. 


PSMB2 

Proteasome  subunit  beta 
type -2  also  known  as  20S 
proteasome  subunit  beta-4 

Forms  a  pivotal  component  for  the  UPS  involved  in  protein 
ubiquitination  and  subsequent  proteolysis  and  degradation 
which  are  important  mechanisms  in  the  regulation  of  the  cell 
cycle,  cell  growth  and  differentiation,  gene  transcription, 
signal  transduction  and  apoptosis.  It  is  also  involved  in 
processing  of  class  I  MHC  peptides 

PSME2 

Proteasome  activator 
complex  subunit  2 

Process  class  I  MHC  peptides 

PSMA5 

Proteasome  subunit  alpha 
type -5  also  known  as  20S 
proteasome  subunit  alpha-5 

Process  class  I  MHC  peptides 

HLADMB 

HLA  class  II 

histocompatibility  antigen, 

DM  beta  chain 

Plays  a  central  role  in  the  peptide  loading  of  MHC  class  II 
molecules  by  helping  to  release  the  CLIP  (class  Il-associated 
invariant  chain  peptide)  molecule  from  the  peptide  binding 
site,  thus  playing  a  major  role  in  MHC  class  II  antigen 
presentation  pathway 

DNMTIA 

DNA  methyltransferase  lA 

Enzyme  catalyzes  the  transfer  of  methyl  groups  to  specific 

CpG  structures  in  DNA  (DNA  methylation).  Considered  to  be 
the  key  maintenance  methyl  transferase  in  mammals, 
predominantly  methylates  hemi  methylated  CpG  di¬ 
nucleotides  in  the  mammalian  genome. 

DNMT3A 

DNA  (cytosine -5)-methyl 
transferase  3A 

DNA  methyl  transferaseresponsible  for  de  novo  DNA 
methylation. 

DNMT3B 

DNA  (cytosine -5-)-methyl 
transferase  3  beta 

DNA  methyl  transferase  responsible  for  de  novo  DNA 
methylation. 

HDACl 

Histone  deacetylase  1 

Class  I  histone  deacetylase,  playing  a  key  role  in  the 
regulation  of  eukaryotic  gene  expression 

HDAC2 

Histone  deacetylase  2 

Class  I  histone  deacetylase,  playing  a  key  role  in 
transcriptional  regulation  and  regulation  of  eukaryotic  gene 
expression 

HDAC4 

Histone  deacetylase  4 

Class  II  histone  deacetylase,  playing  a  key  role  in 
transcriptional  regulation  and  regulation  of  eukaryotic  gene 
expression 
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Primer 

Name 

Forward  Sequence 

Reverse  Sequence 

Product 

Size 

GAPDH" 

TGACAACTTTGGTATYCGTGGAAGG 

AGGCAGGGATGATGTTCTGGAGAG 

134  bp 

ISsrRNA" 

GCTTAATTTGACTCAACACGGGA 

AGCTATCAATCTGTCAATCCTGTC 

69  bp 

PLCEl" 

GCCCAAAGCAAGTGGAAAGG 

TCTTCACCTGGGTTAAACATGC 

700  bp 

IL8 

CAGAGACAGCAGAGCACACA 

GGCAAAACTGCACCTTCACA 

158  bp 

MICE 

CACCCAGGCTGCAGTTCACT 

CGGGAGTCTGAGGTACGAGAA 

88  bp 

PSMB8 

GATCTCCAGAGCTCGCTTTA 

GTTCACCCGTAAGGCACTAA 

200  bp 

CCL5 

CCATGAAGGTCTCCGCGGCAC 

CCTAGCTCATCTCCAAAGAG 

361  bp 

IL6 

GGTACATCCTCGACGGCATCT 

GTGCCTCTTTGCTGCTTTCAC 

88  bp 

IL18 

GCTTGAATCTAAATTATCAGTC 

CAAATTGCATCTTATTATCATG 

335  bp 

DNMTIA 

ACCGCTTCTACTTCCTCGAGGCCTA 

GTTGCAGTCCTCTGTGAACACTGTG 

G 

335  bp 

DNMT3B 

AATGTGAATCCAGCCAGGAAAGGC 

ACTGGATTACACTCCAGGAACCGT 

191  bp 

HDACl 

CCAAGTACCACAGCGATGAC 

TGGACAGTCCTCACCAACG 

110  bp 

HDAC2 

TGAAGGAGAAGGAGGTCGAA 

GGATTTATCTTCTTCCTTAACGTCTG 

124  bp 

DNMT3A 

CGTTGGCATCCACTGTGAATGA 

TTACACACACGCAAAATACTCCTT 

551  bp 

IFNy 

CCAACGCAAAGCAATACATGA 

CCTTTTTCGCTTCCCTGTTTTA 

79  bp 

TNFa 

GGA  GAA  GGG  TGA  CCG  ACT  CA 

CTG  CCG  AGA  CTC  GGG  AA 

70  bp 

ILlp 

GCAAGGGCTTCAGGCAGGCCGCG 

GGTCATTCTCCTGGAAGGTCTGTGG 

GC 

96  bp 

IL4 

ACTTTGAACAGCCTCACAGAG 

TTGGAGGCAGCAAAGATGTC 

74  bp 

IL15 

GTCTTCATTTTGGGCTGTTTCAGT 

CCTCACATTCTTTGCATCCAGATTCT 

316  bp 

TLR2 

GGGTCATCATCAGCCTCTCC 

AGGTCACTGTTGCTAATGTAGGTG 

181  bp 

TLR4 

CAGAGTTGCTTTCAATGGCATC 

AGACTGTAATCAAGAACCTGGAGG 

282  bp 
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PSME2 

GGGAATGAGAAAGTCCTGTCC 

TCAATCTTGGGGATCAGGTG 

113  bp 

IL12 

CCAAGAACTTGCAGCTGAAG 

TGGGTCTATTCCGTTGTGTC 

355  bp 

HDAC4 

GAGAGACTCACCCTTCCCG 

CCGGTCTGCACCAACCAAG 

240  bp 

PSMA5 

AAGCCCATGTTGCTTTTGGG 

GGCGAACGGACATAGGCTAA 

112  bp 

PSMB2 

AGAGGGCAGTGGAACTCCTT 

AGGTTGGCAGATTCAGGATG 

72  bp 

HMGBl 

ACATCCAAAATCTTGATCAGTTA 

AGGACAGACTTTCAAAATGTTT 

122  bp 

HLADMB 

ATGTGAAATCCTTTGGAGTCCCA 

GGCATCTTTACAGAGCAGAGCAT 

145  bp 
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825  “GAPDH  and  ISsrRNA  primers  are  house  keeping  genes.  '’PLCEl  primer  pair  tested  to 

826  amplify  700bp  genomic  region  of  PLCEI  is  used  as  genomic  DNA  control. 
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